A HR-like diagram for galaxies: the M, versus Mqct'^ relation 



A. Feoli 

Dipartimento di Ingegneria, Universita del Sannio, 
Corso Garibaldi n. 107, Palazzo Bosco Lucarelli 
82100 - Benevento, Italy. 

f eoliOunisannio . it 

and 
L. Mancini 

Dipartimento di Fisica "E. R. Caianiello" , Universita di Salerno, 

via S. Allende 
84O8I - Baronissi (SA), Italy. 

ImanciniSphysics . unisa . it 
ABSTRACT 

We show that the relation between the mass of supermassive black holes located 
in the center of the host galaxies and the kinetic energy of random motions of the 
corresponding bulges is a useful tool to study the evolution of galaxies. In the form 
log]^o(-^») = b + mlogiQ{MQa'^ /c^), the best-fitting results for a sample of 64 galaxies 
of various morphological types are the slope m = 0.80 it 0.03 and the normalization 
b = 4.53 lb 0.13. We note that, in analogy with the Hertzsprung- Russell diagram for 
stars, each morphological type of galaxy generally occupies a different area in the M, — 
(Mg<t^)/c^ plane. In particular, we find elliptical galaxies in the upper part of the line 
of best fit, the lenticular galaxies in the middle part, and the late-type galaxies in the 
lower part, the mass of the central black hole giving an estimate of the age, whereas the 
kinetic energy of the stellar bulges is directly connected with the temperature of each 
galactic system. Finally, the values of the linear correlation coefficient, the intrinsic 
scatter, and the obtained by using the M, - Mgct^ relation are better than the 
corresponding ones obtained from the M, — cr or the M, — Mq relation. 



Subject headings: black hole physics - galaxies: general - galaxies: kinematics and 
dynamics - galaxies: statistics 
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Introduction 



Today the fact that many galaxies, of different morphological types, host a supermassive 
black hole (SMBH) at their center has been established on quite solid grounds. The studies of 
the kinematics of galaxies and the combination of multi-band observations have played a ma- 
jor role in this scientific process. At the same time, the idea that the mass of a central SMBH 
is correlated with the evolutionary state of its host galaxy is being consolidated among the sci- 
entific community. In order to qualify this correlation, many relationships have been proposed 
between the mass of the SMBHs and almost all the possible parame t ers of the host galaxy bulge s: 
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199q : iMagorrian et al 



2002; Gebhardt et al 



1998 
200i 



ese &: Merritt 



3 I2QQ0I: 



Gebhardt et al 



Marconi et al 



Marconi &: Hunt 



2001 



Kormendv Sz Richstonel 19951: van der Mare] 19991: 



200E iTrema 



Merritt &: Ferrarese 



2003; iHaring k Rixl 12004 : ICultekin et a 



2001 



ne et al 



200: 



[ lichstone et al 



Laoi 



2001 



galax y light concentration (iGraham et al. 

'"'2OO1I), the X-ray power de nsity spectra 



200 ll ). the dar k matter halo (iFerrarese 2002i'). the radio core length (ICao fc Jianel l2002l '). the ef- 



ner core radius (ILauer et al 
tial (Aller Richstone 



fective radius ( Marconi &: HuntI 2003 ). the Sersic index ( Graham fc Driver 2001, 2007 ). the in- 
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20071 ). the metal abundance ( Kisaka et al I2OO8I ). the core mass deficit (^Kormendv Bender 



2009"), combination of bulge velocity dispersion, effective radius and/ or intensity (Aller Sz Richstonel 



2007i ) . and, very recently, the X-ray l uminosity, the radioluminosity (jGultekin et al.ll2009c]). and ac- 



tive g: alactic nucleus jets (jSokeijl2009l ). An alternative approach has been proposed bv lFeoli Sz Mele 
(|2005l ). who suggested a relation between the black hole (BH) mass and the kinetic energy of ellipti- 
cal galaxies. Then, the se authors extended the study of the relation M, versus Mg(t^ also to lenticu- 
lar and spiral galaxies (jFeoli Sz Meldl2007l ). enlarging their sample to a total of 29 galaxies, and find- 
ing that it has a scatter smaller than the most famous M, — a relation. Here, Mq is the bulge mass 
of the galaxies, where bulge as usual refers to either the hot, spheroidal compo nent (no dark halo 
or di s k contribution) of a spiral/lenticular galaxy or to a full elliptical galaxy (lAller &: Richstone 



20071 : ICattaneo et al.ll2009l ). In particular, for late-type galaxies, Mq is the m ass enclosed within 



a sph ere of radius R fixed from the surface-brightness profile. As an example, lAller Sz Richstone 
(j2007l ) consider that R = W Re is "the best substitute for the bulge mass", with R^ being the 
eff'ective radius of a galaxy. 

Actually, the main problem is that almost all the above-quoted relations are very tight, so 
it is v ery difficult to find, b y studying the scatter of each one of the m, the "most fundamental 
one" (jTremaine et al.l I2OO2I : iNovak et al.l bood : bultekin et al] l2009bl ). Without definitely solv- 
ing this hard problem, the attention of an increasing number of scientists is now focused on the 
M, — a law in order to study the behavior of some peculiar subsets of galaxies. This led to dis- 
cover that the line o f the best fit o f that relationship is different for barred galaxies wi th respect 
to the barless ones (Graham 2008). The sam e occurs for bulges and pseudo -bulges (|H u '2OO8*; 
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a IZhang et al.ll2008l : iGreene k Hoi 12009 ). 



2OO9I). core or coreless (IHul l2008l ) . active or quiescent (jBarth et al.i i2005l : 
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At the same time, from the theoretical point of view, a lot of interesting (analytic al and semian- 



alytical) models were constructed to explain the experimental results (see, for example, 
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)). It is well known that 



a useful method to obtain theoretical predictions, which can be coi npared with th e corre lations de 



rived from experimental data, is based on numerical simulations. iHopkins et al.l (|2007l ) examined 
the origin and the evolution of the correlations between the properties of SMBHs and their host 
galaxies using hydrodynamical simulations of major galaxy mergers, including the effects of gas 
dissipation, cooling, star formation, and BH accretion and feedback. Their simulations suggest the 
existence of a SMBH fundamental plane, analogous to the fundamental plane of elliptical galaxies. 
The best relation that they found (the one with the lowest scatter) is 

logio(M.) = (7.93 ± 0.06) + (0.72 ± 0.12) logio(Afri) + (1.40 ± 0.49) logio(fT20o), (1) 

where M^^ is the galaxy stellar mass in units lO-^-*^ Mq, and (T200 is the bulge velocity dispersion 
in units of 200 km sec~^. These au thors also show the main role played by the kinetic energy 
of random motions first proposed bv iFeoli &: Meld (j2005l ). In particular, they say: "we therefore 
naively expect that the BH mass should scale with M^a'^" , and declare that the correlation between 
the BH mass and the M^kCT^ "is in some sense more basic then the correlation between the BH mass 
and or a" . In other words, their fundamental plane in BH mass can be well represented as a 
"tilted" correlation betw een BH mass a r id the kinetic energy of the random motions in the host 
galaxies (see Figure 10 of I Hopkins et al.l (|2007l )). Another clue is the ratio between the coefficients 
in Equation dH) multiplying the login(Mi* ) and log]^o('^20o) which is very close to 0.5. This is also 
remarked bv lMarulli et al.l (|2008l ). who modeled the cosmological co-evolution of galaxies and their 
central SMBHs within a semianalytic al framework. Their m odel matches well enough the SMBH 
fundamental plane relation derived bv lHopkins et al.l (j2007l ). and their conclusion is identical: the 
SMBH mass does not simply scale with the star formation (stellar mass) or the velocity dispersion 
of the host galaxy. 



The results of iHopkins et al.l (120071 ) and iMarulli et al.l (120081 ) give a strong evidence that 
ga laxy spheroid s and SMBH s do not form and evolve independently and support the approach 
of iFeoli &: Meld (|2005l . 120071 ) . who pointed to the relationship between the masses of the SMBHs 
and the kinetic energy of random motions in their host galaxies. The consequences for the theo- 
retical models of SMBH growth and evolution are non-trivial. 



In the present paper we want to extend the previous analysis of iFeoli &: Meld (120051 . 120071 ) 
(hereafter. Pap er I and Paper II) to a new set of 64 galaxies, almost all extracted by the catalogue of 
GrahamI (j2008l ). The main aim is to probe if the M, — MgO"^ relation is really a helpful instrument 
to study the evolution of the galaxies, that is if it can play the same role as the Hertzsprung- 
Russell (H-R) diagram in the description of the evolution phases of stars. We will see that different 
morphological types of galaxies occupy different positions in the M, — MgO"^ plane, reflecting their 
age and intrinsic features. Finally, we want to confirm that the linear correlation coefficient, the 
intrinsic scatter, and the of our relationship are better than the corresponding values for the 



- 4 - 



M, — 0" or the M, — Mq, law. Our paper is structured as follows. In § 2 we define the samples used 
in our statistics. In § 3 we explain our results and, finally, in § 4 we draw our conclusions. 



2. The samples 



In order to have a homogeneous set of data, we have considered as the main referenc e for the 
masse s of SMBHs and the velocity dispersions of the galaxies the catalogue published by [Graham 
(I2OO8I) . A lternative values (almost compa tible within un c erta inties) can be found for instance in 
Hul ( 20081 ) or in the more recent paper of iGultekin et al. (|2009b ). Our choice involves the values 
of the central velocity dispersions fic, in contrast with our tw o previous papers w here we used the 
effective dispersion velocity a^,. However, as already noted bv lNovak et al.l (j2006l ). the two ways of 
measuring the v eloci t y disp ersion does not generate profound differences. This is also supported 
by the study of IHuI (j2008l ) , who compared the effective dispersion with the central one ar. 



finding that the differences are much smaller than their measurement errors. Also iGultekin et al 



(j2009bl ) compared de to cr^ finding no systematic bias to high or low values. We also remark the 
fact that in our previous papers the data were ex tracted only by sin gle sources: in Paper I all the 
values of the ga. l axies m asses have been taken from lCurir et al.l (119931) , the y elocity dispersions from 



Busarello et al.l (jl992l ) , and the SMBH masses from 



sets of data have been taken from lHaring &: Rixl ( 20041 ) . Here, due to the enlargement of the sample. 



Tremaine et al.l ([2002); in Paper II the three 



pie of data from 



Cappellari et al 



a homogeneous choice is no more poss ible and we are forced to build up a sam 
variou s catalogues - essentially those of lGrahamI (j2008l ) . iHaring &: Rixl (|2004l ) and 
(j200a ) - and single papers. A clear limit of a collection of data of this kind is, of course, related 
to the many different techniques utilized to estimate the masses of the bulges (dynamical or virial 
masses, Schwarzschild models, Jeans equation, etc.; see Appendix [B] for a more comprehensive 
discussion), and of the SMBHs (gas or stellar kinematics, water maser, proper motions, etc.). 

In this paper, we consider two sample s of galaxies. T he first sample (sample A) is composed 
by 49 galaxies included in the table 1 of iGrahamI (|2008l ) . Actually, his catalogue is formed by 
50 galaxies that are considered to have reasonable measurements of their SMB H masses . We 



exclude the galaxy IC2560, s ince a reliable value for its bulge mass is not available (llshihara et al 



2OOII : ISchulz &: Henkelll2003l ). In several cases we would have liked to substitute some data of the 
G raham catalogue with other measures which are less uncertain or simply more recent (like that 
of IGultekin et al.l l2009bl ) , but we did not do so in order to avoid the risk that the tightness of 
our relation might depend on a suitable choice of the data. For example, we have used for the 
mass of SMBH in the Milky Way the val ue cited by Graham even if we know that an update 
value is now available (iGillessen et al.ll2009l ). and we have included in the sample also the elliptical 
NGC221, which we would have liked to exclude from the fit as already done in Papers I and II. 
Starting from the Graham's catalogue we have fixed the total number and the names of galaxies, 
their velocity dispersions, the SMBH masses, and morphological types. In this way, only the galaxy 
mass remains as a free parameter but our choice was anyway restricted by using the data published 
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bv lHaring fc Rixl (|200J) and ICappellari et alJ (120061) 



The second sample (sample B) is composed by t he galaxie s of sa mple A plus other 15 galaxies 
whose parameters have been taken from table 2 of iGrahamI (|2008l ) and from other papers. Of 
course, this enlarged sample does not have the aim to include all the galaxies with a measured BH 
mass or with an upper limit on its value. Sixty-four galaxies are listed in Table [1] and compose 
sample B, whereas only the first 49 galaxies are included in sample A. 

Concerning the errors in the measures, we adopt the same strategy as in Paper II. Following 
Haring &: Rixl (120041 ) . we consider that the error for the bulge mass is 0.18 dex in log^g foi" ^ 
the galaxies, while the relative error on the velocity dispersions is 10%. 



3. Results 



The relation between the mass of the SMBHs and the kinetic energy of random motions of the 
corresponding host galaxies has been presented in Papers I and II in the form 



logio(M.) = 6 + mlogio(MGaVc2). 



(2) 



Thus, this relation can be used to predict the values of M, in other galaxies once we know their 
mass and velocity dispersion. In order to minimize the scatter in the quantity to be predicted, we 
have to perform an ordinary least-squares regression of M, on Mqo^ for the galaxies in Table [U 
of which we already know both the quantities. In Table [21 we compare the fits of our relationship 
for the t wo samples and the corr esponding fits for M, — a and M, — Mq laws. As in Paper II 
(see also lGraham &: Driveii (|2007l )). these fits w ere obtained takin g into account the error bars in 
both variables and using the routine FITEXY (|Press et al.lll992l ) for a relation y = b + mx, by 
minimizing the (see Appendix [A]) . Comparing the results of the three laws, we notice that 
the the intrinsic dispersion eq (i.e., dispersion due to the galaxies themselves rather than to 
measurement errors), and the Pearson linear correlation coefficient r of our relationship are better 
than the other ones (Table [2]). This is also evident by the comparison of Figured^, [lb and Hh, 
where the three relations are reported in log-log plots (we associated a particular marker to each 
galaxy according with its morphological type). It is visually clear that the galaxies (especially the 
spirals) are more spread in the case of the M, — cr or the M, — Mq relations, t han in the M, — Mqu'^ 
law. This result is also supported by the analysis of lGadotti &: KauffmannI (|2009l ) who found that 
elliptical, classical bulges and pseudo-bulges follow different relations between their stellar masses 
and velocity dispersions. 

Comparing the results in Table [2] with the corresponding ones in table 3 of Paper II, we observe 
that, by enlarging the sample, the correlation coefficient of our relationship increases, showing the 
robustness of our idea. Actually, if the exi stence of the correlation is not a novelty, since it has 
already been found bv iFeoli &: Meld (|2005l ). the increase of the correlation coefficients with the 
enlarging of the sample is a result that had not at all been taken for granted. 



As already observed by iNovak et alj (|2006l ) , the question "which relation is better than the 
others?" is extremely sensitive to inaccurate estimates of the measurement errors. So, the result 
that our is better than the M, — a and the M, — Mq laws can be caused by an overestimation 
of the error on the galaxy masses. In order to avoid a similar misleading result, we have checked 
what happens using a standard least squared fitting, assuming that errors in the kinetic energy are 
zero and that errors in the log^o are the same Cy for each galaxy. The results are reported in 
Table [3] and show that the scatter of our relation is better than the M, — a and the M, — Mq laws 
even in this extreme case. Furt hermore, the slope of the line of best fit m = 0.73 it 0.04 is the same, 
inside the errors, as the one of lHopkins et al.l (|2007l ) in Equation ([T]). 



We note also that the slope of the M, — a law depends on the errors and on the fitting 
methods use d, more than t h e oth er rel ations do. Wh i le the values in Table [3] are close to the 
estimates of iTremaine et al.l (l2002l ) and iGultekin et al.l (|2009bl ) (even if they have been obta ined 
with a different sample), the values in Table [2] are closer to the ones obtained by Graham (l2008l^ (he 
found m = 5.22 it 0.40 and b = 8.13 it 0.06) with the same sample but with the lAkritas &: Bershadi 



diaaa) method. 



A surprising result is shown in Figure [21 where we performed a log-log plot of the energy 
stored by the SMBH, E^t = M,c^, as a function of the bulge kinetic energy of random motions, 
both normalized by the rest energy of the Sun, Mqc^. Given the line of best fit (solid line) and a 
sort of border line (dashed line) that divides the diagram in two parts, it is evident that: 



1. almost all the elliptical galaxies (except NGC3377) are in the higher part of the diagram (over 
the dashed line), 

2. the lenticular galaxies are located in the middle-upper part of the diagram, 

3. the barred lenticular galaxies are located in the middle part of the diagram (but under the 
dashed line), 

4. all the spirals are in the middle-lower and in the lower parts of the diagram (under the dashed 
line) , 

5. in the lower part of the diagram we find also two dwarf elliptical galaxies: NGC221 and 
NGC4486A. 



In analogy with the H-R diagram for stars, each morphological type of galaxy occupies a 
different area in the M, - {MGa^)/c^ plane. This effect is in part due to the well known fact that 
M, and Mq (also a, even if with a lot of exceptions) generally increase with the morphological 
type, but it is remarkable and not granted that they simultaneously increase just in the right way 
to produce the effect and a law with a minimal scatter. Compared with the other two relations, 
this clear trend can be lightly recognized also in Figure [Th, but it is not clear in Figure [TJa. For 
example, the ellipticals in Figure [1^ are more separated from the spirals with respect to Figures 
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[TJd or[T}:. This can be quantified calculating the width of the transition area (light red colored), in 
which the elliptical and the spiral galaxies are mixed together, with respect to the entire area (light 
blue colored) occupied by all galaxies (we exclude the dwarf ellipticals). In Figured^, the red zone 
is only the 19% of blue area, whereas it is 33% and 21% in Figure [TJj and[T}:, respectively. 

The general trend observed in Figure [2] is respected in Figure [3] for the galaxies of sample B, 
even if two lenticular galaxies appear in the lower part of the diagram and a spiral galaxy in the 
upper part. The latter is the famous Sombrero gal axy (NGC4594), one of the largest galaxies in 
the nearby Virgo Cluster, classified as a lenticular bv lMagorrian et al.l (|l998l ). It is well known that 
it has a bright nucleus and an unusually large classical bulge, testified by a relatively large number 
of globular clusters. We know that the classical bulges are believed to be generated by mergers 
and are common in early type galaxies but become progressively rare toward later types. They 
share sorne stru ctural, dynamical, and population properties with the lower-luminosity ellipticals 
(lFreemanll2007l ). Actually, NGC4594 is surrounded by a halo of stars, dust, and gas that indicate it 
may actually be described as an elliptical galaxy that contains a more robust interior configuration. 
Therefore, its presence in the upper part of the diagram is not so improper. Later type galaxies 
like the Milky Way mostly have small boxy bulges and are all in the lower part of the diagram. On 
the other han d, both a classical bulge and an inner boxy bulge a re present in NGC224 (Andromeda 
galaxy, M31) ([Athanassoula &: Beatonll2006l : iBeaton et al.ll2007l ). which is located just in the middle 
region. 

Both in Figure [2] and El the elliptical galaxies are all clustered very near the line of best 
fit. Conversely, the galaxies of the other morphological types look slightly more scattered. This 
is particularly true for the lenticular galaxies in the middle-upper part of the diagram. Among 
them, the galaxy NGC4342 is l ocated quite far from the best-fitting line. As already noted by 
Cretton Sz van den BoschI (jl999l ) , NGC4342 is one of the galaxies with the highest SMBH mass to 
bulge mass ratio. The consequent hypothesis that we are in the presence of a galaxy in a particular 
evolu tionary state is also support ed by the presence of both an outer disk and a stellar nuclear 
disk (Ivan den Bosch Sz Jaffd 119971 ). Instead, the position of the lenticular galaxy NGC7457 in the 
lower zone of the diagram is due to its SMBH which is one of the least-massive BHs yet detected 
in the core of a galaxy, roughly the same mass as the BH at the center of our Galaxy. Equally, the 
boxy-bulge lenticular NGC7332 is located in the middle-lower region of Figure [2j Going on, the 
peculiar galaxy NGC5128 (Centaurus A) appears in the middle of the diagram, but still quite far 
from the other lenticulars. The strange morphology of Centaurus A is generally recognized as the 
result of a merger of two smaller galaxies. In this way, it is possible to explain a bulge com prised 
mainly of evolved red stars and a dusty disk, which has been the site of recent star formation (jlsrael 
19981 ). 



Finally, we note the presence of the intermediate-size elliptical NGC3377 in the center of the 
graphic in Figure [21 and three dwarf elliptical galaxies in the lower part. Two of them are NGC4742 
and NGC4486A, both belonging to the Virgo Cluster of Galaxies, whereas the smah NGC221 (M32) 
is a satellite of M31. Continuing the analogies with the H-R diagram, we can look at this area as 
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reserved to the dwarf ellipticals, in the same manner as the region occupied by the white dwarfs 
in a classical color-magnitude diagram. If the three dwarf ellipticals do not really belong to the 
"principal sequence" , we can exclude them from the fit. In this case, the slope and the normalization 
in our relationship for the reduced sample of 61 galaxies are 0.84 it 0.03 and 4.33 ±0.15, respectively 
(the dashed line in Figure [3|). Furthermore, comparing again the Xr = 1-80 and the r = 0.91 of our 
relation with for instance the corresponding Xr — 2-37 and r = 0.83 for the M, — a law, we find an 
increase of the gap between the two relations. 

Also remarkable is the fact that all the barred galaxies (lenticulars and spirals) are located 
only in the lower /middle- lower part of the diagram (under the dashed line of Figure [2|). 

We also study a possible correlation between the activity (Seyfert, Liner, etc.) of each galaxy 
with its position on the diagram, but we did not note any particular trend. 



4. Discussion 

In this paper, we have investigated the relation between the mass of the SMBHs and the 
kinetic energy of the random motion of the corresponding galaxies. This relation has been tested 
on a homogeneous sample of 49 galaxies and then on a more enlarged sample of 64 ones. As 
shown in Table El the statistical analysis confirms the result of our two previous papers, that is the 
proposed relation works well and better than the most common M, — a law or the M, — Mq one. 
Furthermore, the main result that we report consists in the particular positions of the galaxies in 
the M, — (MgO"^)/c^ plane, which resembles the H-R diagram for the stars. Other analogies also 
exist between the two diagrams. The H-R diagram connects the energy radiated (per unit time) 
by the nucleus of a star with its surface temperature. In the same way, our diagram connects a 
property of the inner nucleus of a galaxy, the energy stored by the SMBH, M,c^, with a property of 
the external surface of its bulge, i.e., the kinetic energy of random motions. This energy is related 
just with the temperature of the stellar system. In fact, let us consider a spherically symmetric 
distribution of stars with density p, whose dynamical state is described by a distribution function 
of the form 

^(^) = T:^^^^/^^ (3) 



(2^ct2) 



wher e E = is the binding energy, and ^ is the relative gravitational potential ( Binnev &: Tremaine 



19871 ). Now, it is well known that the structure of a collisionless system of stars, whose density in the 
phase space is given by Equation ([3|), is identical to the structure of an isothermal self-gravitating 
sphere of gas, if we set 

Ma^ = Nk^T, (4) 

with M being the total mass of the system, the number of objects contained in the system, 
k-Q the Boltzmann's constant, and T the temperature of the system. Since a stellar bulge can 
be considered, with a good approximation, similar to a spherically symmetric system, its kinetic 
energy of the random motions Mqu"^ gives an indication of the temperature of the galaxy bulge. 
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On the other hand, since the SMBH at the center of galaxies can only increase its mass, the 
stored energy is directly connected with the initial density of the system and its evolutionary state. 
So, the stored energy of central SMBH will guide the galaxy along the evolutionary process, and, 
in that sense, an accretion of the SMBH, bound up with the flow of time, will imply a migration of 
the galaxy position from the lower-right part of the M, — (Mqct^) diagram to the upper left. This 
migration does not involve a brutal transformation of a spiral galaxy in an elliptical one but, since 
we consider just the mass of the bulge in our relation, we suppose that the spheroidal components 
of the spirals increase their size becoming similar to ellipticals, as it happened for the Sombrero 
galaxy discussed in the previous section. 
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anonymous referee for its suggestions and comments that have helped us to improve the quality of 
this paper. This research was partially supported by FAR fund of the University of Sannio. L.M. 
acknowledges support for this work by MIUR through PRIN 2006 Protocol 2006023491_003, by 
research funds of the Italian Space Agency, by funds of Regione Campania, L.R. n. 5/2002, year 
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A. Appendix 



The formula used in this paper to estimate the maximal errors in the functions F of the 
parameters (a,6, c, ....) is 



AF{a,b, c) 



dF 



da 



Aa + 



OF 



The reduced ) used in Table [H is defined as 

'2 1 



X 



db 



Ab + 



dF 



dc 



Ac. 



(Al) 



N 



N 



N 

2 ^ (Ayi)2 + m2(Axi)2 



(A2) 



for a relation of the form y = b+mx, where is the number of galaxies in the sample. The internal 
scatter eq in Table [2] is calculated for fixed values of m and b imposing that 



mxi 



N 



—y 

-2f^^{AyiY + el + m?{A2 



1. 



(A3) 



Following iTremaine et al.l (|2002l ) , the internal scatter can be computed also replacing the error Ay 
with (Ay2 + Sq)^^^ before the fitting procedure, and then adjusting eo and refitting until the Xr is 
equal to 1. This approach affects the values of the slope and the normalization of the best-fitting 
relations. We reported them in Table H] together with the values of Eq for the galaxies of sample A. 
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Furthermore, the ey used in Tabled is defined as 



N 



4 = 

1=1 



Finally, the Pearson linear correlation coefficient is 

Y17=iixi - x) - iy^ - y) 



(A5) 



B. Appendix 



There are several methods to estimate the mass of a galaxy. Some of them are based on the 
virial theorem, other on the Jeans equation or on self-consistent models. When we are interested in 
the mass of a whole galaxy, then it is possible to use the virial theorem in its scalar form 2T + [7 = 
and, kn owing the rotation ve locity and velocity dispersion, the corresponding virial mass can be 
derived (jBusarello et al.lll99d ). This approach often reduces to calculate the dynamical mass 

Mdyn = (Bl) 

Instead, the Jeans equation is particularly useful to calculate the mass of late type galaxies, 
where we must take into account only the spheroidal component. If one wants to study only the 
part of a galaxy inside a radius R, the system is not completely isolated and t he effect of the matter 
at a radius r > R has to be considered as an unknown external pressure (jChandrasekhaii Il969l ) 
(this term and the rotation velocity are often neglected in the case of elliptical galaxies) . In order to 
overcome this problem, it is convenient to compute the masses of galaxies starting from the Jeans 
equation that describes the equilibrium of a spheroidally symmetric system having eccentricity e 
and an isotropic velocity dispersion tensor: 



1 d 



(r) _ 47rG(l - e^/^ 



dx x'^ p{x) 



(B2) 



pir^dr^^'^ ' " r r (r^ 

where r is the radius in the equatorial (z = 0) plane, p{r) is the (unknown) spatial density, a(r) and 
v(r) are the one-dimensional velocity dispersion and rotation velocity respectively (jBinnev &: Tremaine 
19871 ). The solution of Equation ()B2p can be written in the form p(r) = po x l( r), where Ijr) is the 
lumin osity density. For elliptical galaxies the approach followed for example bv lBusarello &: Longo 
(jl992l ) is to assume that the luminosity distribution corresponds to the spatial deprojection of the 
r^/ "^ law. A simple analytic al approximation for the deprojection of the r^/^ law has been derived 
bv lMellier MathezI (jl983): 

;(r) =r-^exp(-5ri/4), (B3) 

where /? = 0.855 and b = 7.669. Substituting this solution together with V{r) and cr{r) (deprojected 
from the experimental data) in Equation ()B2p . an expression of po as a function of r can be obtained. 
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Computing the value of po for each object in the considered sample of galaxies at 10 different radii, 
the residuals with respect to its mean value < po > turn out to be very small (and will be used 
to estimate the error AM), and show no systematic tren d with the radius, thus s upporting the 
hypothesis that po = constant at least in the inner regions (jBusarello &: Longdll992l ). So, the final 



result for the mass density is p{r) =< po > r^^ exp{—br^^^). Changing the luminosity distribution, 
this method can be applied also to late type ga laxies. A similar approach was followed and well 
explained in their paper by lHaring Rixl ([200J). 



When the galaxy sample is very large and various, it is very difficult that all the masses of 
the galaxies are calculated in the same manner. For t he sample in Table [ll the masses of most of 
the galaxies have been estimated with the pr ocedure of lHaring &: Rixl (|2004l ). a small part with the 
Schwarzschild model (ICappellari et al.ll2006l ). and only in a few cases using the dynamical mass 
formula. 
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Table 1. Sample^' 





lype 






M, 


<5M, 




IVJQ 


JTLtJlcl tJllt/tJo 


CygnusA 


E 


270 


1 


2.5 X 10« 


7.0 X 10* 


1 


1.6 X 10^2 


2 


NGC221 


E2 


72 


1 


2.5 X 10'' 


5.0 X 10"' 


1 


8.0 X 10* 


3 


NGC821 


E6 


200 


1 


8.5 X 10^ 


3.5 X 10' 


1 


1.3 X 10" 


3 


NGC1399 


El 


329 


1 


4.8 X 10* 


7.0 X 10' 


1 


2.32 X 10" 


4 


NGC2974 


E4 


227 


1 


1.7 X 10* 


3.0 X 10' 


1 


1.57 X 10" 


5 


NGC3377 


E5 


139 


1 


8.0 X lO'' 


6.0 X 10^ 


1 


3.08 X 10i° 


5 


NGC3379 


El 


207 


1 


1.4 X 10* 


2.7 X 10* 


1 


6.8 X 10^° 


3 


NGC3608 


E2 


192 


1 


1.9 X 10* 


1.0 X 10* 


1 


9.7 X 10^° 


3 


NGC4261 


E2 


309 


1 


5.2 X 10* 


1.1 X 10* 


1 


3.6 X 10" 


3 


NGC4291 


E2 


285 


1 


3.1 X 10* 


2.3 X 10* 


1 


1.3 X 10" 


3 


NGC4374 


El 


281 


1 


4.64 X 10* 


3.46 X 10* 


1 


3.6 X 10" 


3 


NGC4473 


E5 


179 


1 


1.1 X 10* 


8.0 X lO'' 


1 


9.2 X 10^° 


3 


NGC4486 


EO 


332 


1 


3.4 X 10^ 


1.0 X 10^ 


1 


6.0 X 10" 


3 


NGC4486A 


E2 


110 


1 


1.3 X 10'^ 


8.0 X 10^ 


1 


4.06 X 10^ 


6 


NGC4621 


E5 


225 


1 


4.0 X 10* 


6.0 X 10' 


1 


1.88 X 10" 


5 


NGC4649 


El 


335 


1 


2.0 X 10'-^ 


6.0 X 10* 


1 


4.9 X 10" 


3 


NGC4697 


E4 


174 


1 


1.7 X 10* 


2.0 X 10' 


1 


1.1 X 10" 


3 


NGC5077 


E3 


255 


1 


7.4 X 10* 


4.7 X 10* 


1 


2.1 X 10" 


7 


NGC5813 


El 


239 


1 


7.0 X 10* 


1.1 X 10* 


1 


5.05 X 10" 


5 


NGC5845 


E3 


233 


1 


2.4 X 10* 


1.4 X 10* 


1 


3.7 X IQi" 


3 


NGC5846 


EO 


237 


1 


1.1 X 10^ 


2.0 X 10* 


1 


6.36 X 10" 


5 


NGC6251 


E2 


311 


1 


5.9 X 10* 


2.0 X 10* 


1 


5.6 X 10" 


3 


NGC7052 


E4 


277 


1 


3.7 X 10* 


2.6 X 10* 


1 


2.9 X 10" 


3 


NGC3115 


SO 


252 


1 


9.1 X 10* 


1.03 X lO'' 


1 


1.2 X 10" 


3 


NGC3245 


SO 


210 


1 


2.1 X 10* 


5.0 X lO'' 


1 


6.8 X 10^0 


3 


NGC3414 


SO 


237 


1 


2.5 X 10* 


4.0 X 10^ 


1 


1.7 X 10" 


5 


NGC3998 


so 


305 


1 


2.2 X 10* 


2.0 X 10* 


1 


5.5 X 10^° 


8 


NGC4342 


SO 


253 


1 


3.3 X 10* 


1.9 X 10* 


1 


1.2 X IQi" 


3 


NGC4459 


SO 


178 


1 


7.0 X 10^ 


1.3 X lO' 


1 


7.86 X 10^° 


5 


NGC4552 


so 


252 


1 


4.8 X 10* 


8.0 X lO'' 


1 


1.87 X 10" 


5 


NGC4564 


so 


157 


1 


5.6 X lO'' 


3.0 X 10« 


1 


4.4 X IQi" 


3 


NGC5128 


so 


120 


1 


4.9 X lO'' 


1.8 X lO'' 


1 


2.16 X 10i° 


9 


NGC5252 


so 


190 


1 


1.06 X 10^ 


1.63 X 10^ 


1 


2.4 X 10" 


2 


NGC1023 


SBC 


204 


1 


4.4 X 10'' 


5.0 X 10"^ 


1 


6.9 X IQi" 


3 


NGC2778 


SBO 


162 


1 


1.4 X 10'^ 


9.0 X 10'3 


1 


1.06 X lOio 


10 


NGC2787 


SBO 


210 


1 


4.1 X 10'^ 


5.0 X lO'' 


1 


2.9 X IQi" 


11 


NGC3384 


SBO 


148 


1 


1.6 X lO'^ 


2.0 X 10'3 


1 


2.0 X 10^° 


3 


NGC4596 


SBO 


149 


1 


7.9 X 10'^ 


3.8 X lO'' 


1 


2.6 X IQio 


2 


Circinus 


s 


75 


1 


1.1 X lO*' 


2.0 X lO'' 


1 


3.0 X 10'' 


12 


NGC224 


SA 


170 


1 


1.4 X 10* 


9.0 X 10' 


1 


4.4 X 10"' 


13 


NGC3031 


SA 


162 


1 


7.6 X 10' 


2.2 X 10' 


1 


1.0 X IQi" 


14 


MW 


SB 


100 


1 


3.7 X 10** 


2.0 X 10^ 


1 


1.1 X lO^" 


3 


NGC1300 


SB 


229 


1 


7.3 X 10' 


6.9 X 10' 


1 


2.14 X 10i° 


15 


NGC3079 


SB 


146 


1 


2.4 X 10"^ 


2.4 X 10'^ 


1 


1.7 X 10^ 


16 


NGC3227 


SAB 


133 


1 


1.4 X lO'' 


1.0 X lO'' 


1 


2.95 X 10^ 


17 
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Table 1 — Continued 



Galaxy 


iype 


0"c 


References 


/l/f 




References 


Mq 


References 






(km/s) 




{Mq) 


[Mq] 








NGC4151 


SAB 


156 


1 


6.5 X lO'^ 


7.0 X 10® 


1 


1.09 X 10" 


18 


NGC4258 


SAB 


134 


1 


3.9 X lO'^ 


1.0 X 10® 


1 


1.1 X 10^^ 


2 


NGC4945 


SB 


100 


1 


1.4 X 10'^ 


1.4 X 10® 


1 


3.0 X 10^ 


12 


NGC7582 


SB 


156 


1 


5.5 X lO'' 


2.6 X lO'^ 


1 


1.31 X 10" 


19 


IC1459 


E3 


306 


1 


1.5 X 10^ 


1.0 X 10^ 


2 


6.6 X 10" 


2 


IC4296 


E 


336 


1 


1.3 X lO'' 


4.0 X 10* 


1 


1.56 X 10^2 


20 


NGC3607<= 


E 


229 


21 


1.2 X 10** 


4.0 X lO'^ 


21 


2.70 X 10" 


21 


NGC4486B 


EO 


169 


1 


6.0 X 10* 


3.0 X 10* 


1 


1.22 X 10" 


22 


NGC4742 


E4 


109 


1 


1.4 X lO'' 


5.0 X 10® 


1 


6.2 X 10^ 


3 


NGC5576<= 


E3 


183 


21 


1.8 X 10® 


4.0 X 10^ 


21 


1.47 X 10" 


21 


NGC3585<= 


SO 


213 


21 


3.4 X 10* 


1.5 X 10* 


21 


1.85 X 10" 


21 


NGC4026<= 


SO 


180 


21 


2.1 X 10* 


7.0 X 10^ 


21 


5.17 X IQio 


21 


NGC7332 


SO 


135 


1 


1.3 X 10^ 


6.0 X 10® 


1 


1.5 X 10^° 


3 


NGC7457 


SO 


69 


1 


3.5 X lO® 


1.4 X 10® 


1 


7.0 X 10^ 


3 


NGC4203 


SBO 


124 


11 


5.2 X lO'' 


1.0 X 10® 


11 


1.5 X 10^° 


11 


NGC1068 


SA 


151 


1 


8.4 X 10^ 


3.0 X 10^ 


1 


1.5 X 10^° 


23 


NGC2748 


SA 


79 


24 


4.4 X 10^ 


3.6 X lO'' 


15 


1.69 X IQi" 


15 


NGC4594 


SA 


240 


1 


1.0 X lO'' 


1.0 X lO'' 


3 


2.7 X 10" 


3 


NGC7469 


SAB 


152 


25 


1.22 X 10^ 


1.40 X 10® 


25 


4.5 X 10^ 


26 
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Note. — ^Adopting the same strategy as in Paper II, we consider that the error for the bulge mass is 0.18 dex in 
login ^^G for '^ll the galaxies, while the relative error on the velocity dispersions is 10%. ''Galaxy types are taken from 
l&aham (2008) with the following exceptions: IC4296, NGC221, NGC821, NGC1399, NGC2974, NGC3379, NGC3607, 
NGC4374, NGC4486A, NGC4486B, NGC4621, NGC5077, NGC5813, NGC5846, NGC6251, which are taken from the 
NASA/IPAC Extragalactic database. "^Effective dispersion Ce is given rather than the central ctc- 
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Table 2. Fitting parameters of the M, — Mgct^, M, — a and M, — Mq relations, both for sample 

A and B 



Relation 


Sample 


N 


m ± Am 


b±Ab 


X? 


£0 


r 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


M. - Mqct'^ 


A 


49 


0.80 ±0.03 


4.49 ±0.15 


1.74 


0.19 


0.92 


M, - cr 


A 


49 


5.06 ± 0.25 


8.18 ±0.04 


1.85 


0.25 


0.87 


M. - Mq 


A 


49 


1.15 ± 0.05 


-4.35 ± 0.57 


2.10 


0.27 


0.89 


M, - Mgo-2 


B 


64 


0.80 ±0.03 


4.53 ±0.13 


1.92 


0.21 


0.92 


M. ~ a 


B 


64 


5.00 ±0.21 


8.20 ±0.04 


2.34 


0.32 


0.85 


M. - Mq 


B 


64 


1.13 ± 0.05 


-4.19 ± 0.49 


2.77 


0.26 


0.90 



Note. — We report the used relation in Column 1, the sample in Column 2, and 
the corresponding number of galaxies in Column 3. By using the routine FITEXY, 
we find the best fit of the relationship y = b + mx. The results for m and b are in 
Columns 4 and 5 and the corresponding reduced x? = /{N — 2) in Column 6. The 
internal (intrinsic) scatter (Column 7) is estimated as that which yields a x? = 1 with 
respect to the given best- fitting relation (see also Appendix [XJ . Finally, the linear 
correlation coefficient is shown in Column 8. 



Table 3. Fitting parameters of the M, — Mqu'^, M, — a and M, — Mq relations, both for sample 
A and B, without considering errors in the kinetic energy 



Relation 


Sample 


N 


m ± Am 


b±Ab 


4 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


M. - Mgo-^ 


A 


49 


0.74 ± 0.04 


4.80 ±0.20 


0.10 


M, - a 


A 


49 


4.46 ± 0.36 


8.13 ±0.06 


0.16 


M, - Mg 


A 


49 


0.98 ±0.07 


-2.56 ±0.79 


0.15 


M, - Mgo-2 


B 


64 


0.73 ± 0.04 


4.88 ±0.18 


0.11 


M. - cr 


B 


64 


4.12 ±0.32 


8.17 ±0.06 


0.19 


M, - Mg 


B 


64 


0.99 ±0.06 


-2.60 ±0.67 


0.14 



Note. — We report the used relation in Column 1, the sample in 
Column 2, and the corresponding number of galaxies in Column 3. By 
using a standard least squared fitting and assuming that errors in the 
galaxy mass and velocity dispersion are zero and that errors in the 
logj^Q Mm are the same ey for each galaxy, we find the best fit of the 
relationship y = b + mx. The results for m and b are in Columns 4 
and 5, respectively. The values of the ey are in Column 6. 
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Table 4. Fitting parameters of the M, — Mgct^, M, — a and M, — Mq relations, for sample A, 
with a different method to calculate the internal scatter 



Relation 


Sample 


N 


£0 


m ± Am 


6± A6 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


M. - Mgo-^ 


A 


49 


0.110 


0.78 ± 0.04 


4.61 ± 0.19 


M. - a 


A 


49 


0.230 


4.25 ± 0.35 


8.18 ± 0.06 


M, - Mq 


A 


49 


0.165 


1.07 ±0.07 


-3.55 ±0.78 



Note. — We report the used relation in Column 1, the sample in 
Column 2, and the corresponding number of galaxies in Column 3. The 
values of the internal scatter in Column 4 are calculated with the proce- 
dure explained in Appendix |X] The corresponding values for the slope 
and the normalization are reported in Columns 5 and 6, respectively. 
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Fig. 1. — Best-fitting (a) M, - Mca'^, (b) M, - a, and (c) M, - Mq relations for the elliptical 
galaxies (red ellipses), lenticular galaxies (green circles), barred lenticular galaxies (green barred 
circles), spiral galaxy (pink spirals), barred spiral galaxies (blue barred spirals), and dwarf elliptical 
ffa.la.xies ('ora.ne'e rmind ellinses^ nf Ra.mnle A. The lierht-red cnlnred area, renresents the tra.nsitinn 
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Fig. 2. — For the M, — Mqu^ relation, the hne of best fit for the galaxies of sample A is represented 
by a solid line, while the dashed line separates the early type galaxies from the late type ones. The 
markers are the same of Figure [TJ 
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Fig. 3. — For the M, — Mqct^ relation, the hne of best fit for the galaxies of sample B is represented 
by a solid line, whereas the one obtained without considering the three dwarf ellipticals by a dashed 
line. The markers are the same of Figure [H 



